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ULTRAVIOLET EXTINCTION AT HIGH GALACTIC LATITUDES 
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ABSTRACT 

In order to study the properties and effects of high Galactic latitude dust we present an analysis of 
373,303 galaxies selected from the Galaxy Evolution Explorer (GALEX) All-Sky Survey and Wide- 
Field Infrared Explorer (WISE) All- Sky Data Release. By examining the variation in aggregate 
ultraviolet colors and number density of these galaxies we measure the extinction curve at high latitude. 
We additionally consider a population of spectroscopically selected galaxies from the Sloan Digital 
Sky Survey (SDSS) to measure extinction in the optical. We find that dust at high latitude is neither 
quantitatively nor qualitatively consistent with standard reddening laws. Extinction in the FUV and 
NUV is ~10% and ~35% higher than expected, with significant variation across the sky. We find 
that no single Ry parameter fits both the optical and ultraviolet extinction at high latitude, and that 
while both show detectable variation across the sky, these variations are not related. We propose that 
the overall trends we detect likely stem from an increase in very small silicate grains in the ISM. 

Subject headings: astronomical databases: atlases, ISM: dust, extinction. Galaxy: local interstellar 
matter, galaxies: photometry 



1. INTRODUCTION 

The interstellar medium (ISM) is suffused with astro- 
physical dust, a fine particulate material ranging from 
plentiful large molecules up to rare micron-sized grains 
([Draind |2QQ3[ ). This dust is important both for the 
physics of the formation of planets and stars in the mod- 
ern universe, as well as mediating the thermal charac- 
teristics of the ISM. It also plays a crucial role in our 
observations of everything beyond the solar system, af- 
fecting light through scattering, absorption, and emission 
from radio wavelengths through X-rays. Understanding 
the distribution, structure, and properties of dust is thus 
a central problem in modern astronomy. 

Extinction by dust at low opacity {Ay < 1) is es- 
pecially important for the study of objects beyond the 
Milky Way. We are now observing the cosmos with 
extremely high photometric precision in regions of the 
sky with low extinction. We have a particular inter- 
est in the effects of high-latitude dust in the context 
of large-area surveys, as they cover many disparate re- 
gions of the sky but aim to compile consistent samples 
of objects. Specifically, errors in dust extinction can 
have disastrous effects on the determination of cosmo- 
logical parameter s (e. g. Huterer, Cunha, & Fang, 2013, 
[Menard et al. 2QlQaj). Secondarily, the study of diffuse 
dusty gas outside of galaxies is being explored for the 
first time, demanding a better understanding of how the 
observed extinction compares to t he actual quantities of 
dust in low density enviro nments (jMenard et al.l l2QlQbl 
iMenard fc Fukugital[2QT2l ) . More detailed studies of dust 
in these low-density regimes are crucial for these studies. 

A great advance was made by 

iSchlegel et all (J1998I : SFD) who provided a Galac- 
tic dust map, using all-sky surveys in the far infrared 
by IRAS and DIRBE. Under the assumption of a 
single, consistent grain population, SFD employed the 
optically thin 100 micron flux from thermally emitting 
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dust grains to estimate the column of dust along 
the line of sight, correcting for variation in observed 
temperature in the dust. Recently, there has a been 
a burst of activity mining Sloan Digital Sky Survey 
(jYork et al.l 12000) to make data sets of "standard 
crayons" (objects whose observed aggregate color vari- 
ations indicate foreground reddening) to study dust at 
high Galactic latitude and test the accuracy of SFD. 
Schlafiv et al. ('2010') used the blue tip of the stellar locus. 



Peek fc Gr aves (2010; PGIO) used spectrosco pic observa- 



tions of quiescent non-starforming galaxies, iJones et"al] 
(120111) used spectroscopic ob servations of M-dwarfs, and 
iSchlafiv fc Finkbeinerl (|2011l : SFll) used spectroscopic 
observations of a broad range of stars. All of these 
explorations found that SFD is largely accurate in the 
optical, but that overall errors and spatial variability of 
the reddening are detectable. In particular, SFll found 
that SFD overes timates reddening by 14% and that the 
reddening l aw of lFitzpatric3 (119991; F99) is favo red over 
lO^Donnelll ([1991 O^D94) or lGardelh d~an ([1989; CCM). 
We extend these explorations of high-latitude extinc- 
tion into the ultraviolet using GALEX observations 
of galaxies as our standard crayons. The ultravio- 
let is an inherently more sensitive waveband for the 
study of dust, as it is typically ~ 3 times more extin- 
guished than the optical. Furthermore, the GALEX 
near-UV channel (NUV) covers the 2175 A bump, 
possibly a probe polycycli c aromatic hydrocarbons 
(PAHs; ISteglich et aP [20Tol ). The GALEX far-UV 
channel (FUV) probes the far-UV rise of the extinc- 
tion curve, associated with v ery small silicate grains 
(jWeingartner fc Draind l2001at WDOl). Small grains 
dominate the heating of the diffuse ISM (jWolfire et al.l 
l2003f ) and may act as tracers of variation in grain size 
distribution as a whole. Thus these grains are important 
both for the physics of the ISM, and better understand- 
ing the variation in grain size. CCM, 0'D94, and F99 all 
report reddening laws in which the average UV extinc- 
tion can be predicted from the overall dust column and a 
single parameter, Ry = Ay/E {B — V). These laws are 



derived by comparing stars extinguished by Ay ^ 1 to 
relatively unextinguished stars, and examining the rela- 
tive colors. It is standard practice to extrapolate these 
extinction laws to the vast areas of the sky of much lower 
extinction in which they have not been adequately tested. 
We rectify that problem here by measuring extinction at 
high latitudes directly. 

In ^we describe the observations and selection criteria 
that go into making our UV standard crayon sample of 
galaxies. In ^we lay out how we grid the galaxy samples 
and how we compare these grids to Monte Carlo models. 
In ^ we report the results of the analysis and provide 
sky maps of color, reddening and extinction. We discuss 
the implications of these results in ^ and conclude in 
^ We include an Appendix that details some tests of 
systematic errors in our data and SFD. 

2. OBSERVATIONS AND DATA 

In addition to the published SFD reddening maps and 
the standard crayon catalog of quiescent galaxies from 
PGIO, we use the full sky photometric catalogs from the 
GALEX satellite and the WISE satellite. 

2.1. GALEX 

GALEX data were ob tained as part o f the GALEX 
All-Sky Survey (AIS) ([iMartin et al.(l2QQ5D. processed us- 
ing the standard GALEX pipeline (|Morrissev et al.ll2QQ7l ) 
and released in catalog form as part of GCAT (Seibert 
et al. , in prep)^ based on the GALEX GR6 data release. 
We selected ah GALEX objects in the GCAT ASC table 
with NUV < 20, a total of 7,577,668 objects with median 
exposure time in FUV and NUV of 144s. The objects 
must be observed in GALEX FUV ^ although not nec- 
essarily detected. Because GALEX all-sky observations 
included in GR6 avoided the Galactic plane, these data 
only include a small number of objects (< 5%) at low 
Galactic latitude \b\ < 10. Standard GCAT aperture 
magnitudes and errors were used for FUV and NUV. 

Although we are using standard GALEX products we 
highlight one potential source of systematic error that 
could influence the results reported below. We are most 
concerned with any systematic in FUV — NUV^ which 
could result from a z ero-point offset in eith er band. Us- 
ing WD standards, (|Morrissev et al.l [20071 ) report ran- 
dom magnitude errors ±0.05, ±0.03 for FUV, NUV re- 
spectively, and drifts ^ 0.015 yr~^ and systematic correc- 
tions between releases ~ 0.01 — 0.02. This systematic is 
smaller than the effects we describe below by factors of a 
few, suggesting that future re-calibrations of the GALEX 
zero-point could have a quantitative but probably not a 
qualitative impact on the results we present here. 

2.2. WISE 

We obtain Wid e-Field Infrared Explorer (WISE; 
iWright et ani2010[ ) catalogs from the All-Sky Data Re- 
lease (12 March 2012) Source Catalog archived at IRSA 
(NASA/IPAC Infrared Science Archive). This cata- 
log contains photometry for objects detected as a point 
source (5cr) in any of the four WISE bands (WISE 3.4, 
4.6, 12, 22 micron bands; Wl, W2, W3 and W4). We ex- 
tract WISE photometric measurements for objects with 
positions within 4 arcseconds of the positions of GALEX 
A/" t/V- detected objects described above. We require the 



galaxies to be detected in Wl, W2, W3 in areas of sky 
observed > 8 times (> 60 s), the nominal WISE All-Sky 
Survey depth. The WISE PSF is comparable to that of 
GALEX, - 6 arcsecond (FWHM). Ah WISE photometry 
discussed in this paper are from WISE pipeline profile-fit 
fluxes. We assume that all of our objects are unresolved 
by WISE, and have verified that our results are insensi- 
tive to this assumption. 

2.3. Galaxy Samples 

We use two galaxy samples to act as standard crayons, 
which allow us to probe high-latitude Galactic extinc- 
tion along the line of sight in both the optical and the 
ultraviolet. 

The first galaxy sample we develop is taken from 
the WISE-selected sample of objects that appear in the 
GALEX AIS photometric data set (Figure |T]). We then 
make color and magnitude cuts based on the WISE and 
GALEX data to select for galaxies (see Figure [1]). Firstly, 
we exclude objects with NUV > 20, as their color distri- 
bution is too broadened by noise to be useful at the depth 
of the AIS. A magnitude cut at W3 < 11.5 finds only ob- 
jects with SNR > 8 in the W3 band. A color cut at W3 - 
W2 > 2 removes some stars from the sample and a color- 
color cut of (W2 - W3-3.15)/2 > (Wl - W2 - 0.76) 
eliminates quasars (IWu et al.l 12011 ). The final G ALEX- 
WISE galaxy sample includes 373,303 galaxies. We note 
that none of our IR magnitude or color cuts are de- 
pendent upon foreground extinction, as extinction in 
the mid-IR at the latitudes where AIS data are taken 
is negligible. As a check, we compare the ultraviolet 
color of these galaxies with GALEX AIS objects with 
NUV < 20 matched to the full spectroscopic SDSS MGS. 
These SDSS galaxies show an almost identical distribu- 
tion in FUV — NUV color. When visually examining 
SDSS photometry of the subset of GALEX- WISE ob- 
jects that reside in the SDSS photometric area we find 
only slightly extended objects consistent with galaxies. 
It is simple to take the complementary WISE color-color 
cut ((W2 - W3 - 3.15) /2 < (Wl - W2 - 0.76)) to gen- 
erate a sample of 95,846 quasars. While this sample is 
smaller and somewhat broader in FUV — NUV color, 
it serves as a sanity check against the GALEX- WISE 
galaxy sample. 

The second data set is essentially identical to the 
151,637 galaxies described in PGIO. These galax- 
ies are spectrosc opically selected from the SDSS Main 
Galaxy Sample (iStrauss et al.l [2002f ) Data Release 7 
(jAdelman-McCarthv et al.l l2008l ) to have neither de- 
tectable H a nor [Oil], and thus be non-starforming, 
using the v alues from the NYU Value Added Catalog 
(jBlanton et al. 2005). Their colors are then corrected 
for the color-magnitude relation and color-density rela- 
tions, as functions of redshift. As a result, their Galactic 
reddening-corrected colors have a extremely narrow dis- 
tribution, with typical scatter of only 25 millimagnitudes 
in ^f — r color. See PGIO for the details of this method. 
The only way the data set we use differs from the PGIO 
galaxies is that in this work we examine the colors of the 
galaxies uncorrected for Galactic extinction. 

3. METHODS &; ANALYSIS 

Our goal is to investigate the dependency of ultravi- 
olet and optical reddening on the characteristics of the 
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Figure 1. Selection criteria for the G ALEX- WISE galaxy sam- 
ple. Green points represent the entire sample of objects matched 
between GALEX AIS and the WISE All-Sky survey. Red points 
represent objects with W3 < 11.5 (red line cut shown in lower 
left panel), blue points further meet the W3 - W2 > 2 and 
(W2 - W3 - 3.15) /2 > (Wl - W2 - 0.76) criteria (blue hue cut 
shown in upper left panel). The lower right shows the similarity 
between the resultant galaxy sample (blue) and the SDSS MGS 
sample (black), in terms of FUV — NUV color. Only small, ran- 
domly selected sub-samples are shown to reduce confusion, and 
match the sample sizes in the bottom-right panel. 

intervening gas, specifically the expected reddening along 
the line of sight determined from the far IR emission of 
the grains, as measured by SFD. This grain column is 
reported by SFD98 as an optical reddening, E{B — V); 
SFll demonstrated that the actual value of the redden- 
ing in these bands may be somewhat smaller than was 
originally proposed in SFD. To avoid confusion, through- 
out this work we will refer to the original value reported 
by SFD as £^sfd, which readers might find easier to inter- 
pret as an extinction, £^sfd = ^v,sfd/3.1. To compare 
our observed extinction and reddening to extant redden- 
ing laws we use a modeled spectral energy distribution 
(SED) for a collection of 4000 galaxies that meet our se- 
lection criteria described in ^2.3[ which h as a very modest 
effect as compared to a simply fiat SED (|da Cunha et alJ 
[2003). 

In addition to simply determining the aggregate ex- 
tinction and reddening parameters in the UV channels 
for the whole of the GALEX- WISE sky (ED) we wish 
to explore the variation and structure of ultraviolet ex- 
tincting dust. We develop the necessary methods below. 

3.1. Gridding 

Were our galaxies to be all identical copies, with zero 
intrinsic variation in color, and were our data to be noise- 
free, we could examine the relationship between color and 
^SFD by simply showing a scatter plot. Nature is, of 
course, neither so simple nor so boring, and the intrinsic 
variation in both of our galaxy samples makes finding 
departures from the expected reddening along lines of 
sight to individual galaxies impossible. To explore the 
structure of ultraviolet extinction we bin the data by 
position on the sky, making maps with pixels of 16 square 
degrees. This is a sensible approach in that it allows us 
to make maps of the color variation and provides a data 
set of points with well-determined ^sfd, optical color, 
and UV color values. We use this pixel area to strike a 
balance between sensitivity and resolution. 

To construct these maps we bin the data on a two 



Zenithal Equal Area grids, one for each Galactic hemi- 
sphere. An equal area projection is useful as it simplifies 
the significance variation in the data from pixel to pixel, 
and a projection center towards the zenith (and nadir, in 
the case of the Southern Galactic cap), reduces distor- 
tion in a map with little data in the Galactic plane. We 
note that it is more common to use Hammer- Aitoff pro- 
jections for Galactic maps, but these tend to maximally 
distort at the poles. 

To determine the galaxy color assigned to that position 
in the map, we take the median of all galaxy colors within 
a pixel, Ca-b' We compute errors for each pixel in two 
ways and use the maximum of the two to be conservative. 
First, we compute a standard bootstrap error for the 
median color in each pixel, a^, resampling the galaxies 
with replacement. Second, we determine a typical width 
for the overall distribution by fitting Gaussian to the 
color histogram, a\ and generate errors for each pixel 
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This method avoids problems with low-number statistics 
in bootstrap errors. There are two differences between 
the PGIO maps and the ones we make in this work. The 
pixels in our maps are statistically independent, and have 
no weighting function. The PGIO maps use a smooth- 
ing kernel to aid interpolation, and thus report a PSF- 
weighted median color and have covariant errors. This 
method can be applied to any quantity evaluated at the 
position of the galaxies. A key point is that we evalu- 
ate SFD maps at specific locations (the galaxy positions) 
and then rebuild maps from those locations. We choose 
to follow this procedure because it allows us to preserve 
the weighting in each pixel defined by the distribution of 
galaxies that contribute to the color map. If we were to 
use a map that was not based to the distribution of galax- 
ies, we would be sampling observed reddening and £^sfd 
from slightly different areas of sky, which could lead to 
systematic errors, especially given the natural bias of a 
magnitude-limited sample of galaxies toward regions of 
lower extinction. We refer to this galaxy weighted me- 
dian as £^SFD- 

Implicit in our analysis is that galaxies act as identi- 
cal objects; standard crayons in the terminology of PGIO. 
While the selected galaxies vary significantly in color and 
many other properties, we are only concerned whether 
they intrinsically vary as a function of their position on 
the sky. In PGIO we found there was a slight trend that 
the selected galaxies were somewhat redder in denser en- 
vironments, and corrected for this trend. Here we exam- 
ine a similar metric, S5, defined as the surface density of 
galaxies angularly closer than the 5th nearest neighbor- 
ing galaxy. We evaluate this metric over low extinction 
regions of the sky, and find no detectable correlation be- 
tween this measure and UV galaxy color. 

3.2. Modeling 

Maps of galaxy color and associated £^sfd allow us to 
make some inferences about the variation in optical and 
ultraviolet reddening across the high-latitude sky (see 
^4.2p . but to learn about Rnuv = Anuv/Esfb and 
Rfuv = ^Ff/y/^SFD directly, we must take into ac- 
count both the color and the number density of galaxies 



on the sky. Note that here we are defining the Ruv^ with 
respect to ^sfd, rather than Eb-uv- To do this we con- 
struct a suite of simple but complete Monte Carlo simu- 
lations of the UV sky. For a choice of Rpuv and Rnuv 
we build a full mock distribution of galaxies on the sky, 
which we then grid as above in §3.11 First, we construct a 
subsample our 373,303 galaxies on the sky that have very 
low reddening, 0.02 < £^sfd < 0.03, which we refer to as 
the comparison sample (CS; see Appendix for discussion 
of errors in SFD at £;sfd < 0.02). We use the CS, along 
with the intersection of the GALEX AIS exposure map 
and the WISE 8-exposure map, to determine the density 
of galaxies on the sky at low extinction, rig. We then 
populate each map pixel area with galaxies drawn ran- 
domly from this distribution, matching the mean density 
observed at low extinction, as well as an empirically mo- 
tivated scatter as described in Equation [2l below. The 
large scale structure (LSS) of the universe causes signif- 
icant variation in galaxy number counts than would be 
expected from simple Poisson noise. A full accounting 
of the 2-point (and N-point) correlation function of the 
galaxies in this sample is beyond the scope of this work, 
however, by examining the low-extinction region of the 
sky from which we drew the CS, we can determine the 
variation in galaxy counts as a function of the fraction 
of the sky observed in a given pixel. We find that the 
standard deviation in the number of galaxies in a pixel 
of area ^ is well fit by the simple function 
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While this method does not take into account any co- 
variation in galaxy count of neighboring pixels, none of 
our results seem to depend on modeling this covariation. 

Model galaxies are de-extinguished according to their 
original observed £^sfd and re-extinguished according to 
the £^SFD observed in their placed position and the cho- 
sen Rfuv a nd R nuv- The data are then gridded as 
described in ^3.11 We also count the number of galaxies 
in each pixel Nnuv^ and the number of galaxies in each 
pixel that have FUV < 20, A^fuv, as a direct constraints 
on Rnuv and Rfuv^ respectively. We note that the ex- 
act way in which it is appropriate to de-extinguish the 
subsample galaxies is unclear. Is the 0.02 < £^sfd < 0.03 
extinction reliable? Should we assume that i^FUV and 
Rnuv are the canonical values or consistent across the 
sky? Is there a significant enough issue with the £^sfd 
zero-point such that de-extincting the galaxies at low ex- 
tinction at all is too unreliable? While these effects do 
quantitatively effect our results weakly at very low ex- 
tinction, they have no qualitative effect on our conclu- 
sions. 

The Monte Carlo model is then run 1000 times for each 
pair in a broad range of Rfuv and Rnuv- We compare 
the observed values of Cfuv-nuVi NnuVi and Nfuv 
for each pixel to the cumulative distribution functions of 
these same statistics in a given Monte Carlo run. Tak- 
ing into account covariance between these three statis- 
tics, we determine \/Xred ^^^ each pixel for each model. 
This is equivalent to finding for each pixel where in the 
multivariate tolerance region of the model the observed 
Cfuv-nuVi Nnuv^ and Nfuv he. Thus we determine 
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Figure 2. The GALEX- WISE galaxy distribution as a function of 
^SFD- The top panel shows the cumulative distribution function 
(CDF) of the GALEX- WISE galaxy sample and the CDF of the 
whole sky, each as a function of £^SFD- The middle panel shows 
the number density of the galaxies and the comparison sample, CS. 
The bottom panel shows the measured Rnuv from the number 
density ratio to the CS and a second-order polynomial fit and 95% 
confidence contours. Also shown is the expectation from 0'D94 
and F99 reddening laws. 

what range of Rfuv and Rnuv provides an acceptable 
fit for each pixel in our maps. 

4. RESULTS 
4.1. Overall Exctinction Parameters 

To determine typical values for Rfuv and Rnuv over 
the whole of the GALEX- WISE sky, we simply bin galaxy 
data by £^sfd , and measure colors and number counts of 
galaxies within those bins. We note that the area from 
which the GALEX- WISE galaxies were drawn represents 
nearly all of the sky at £^sfd — 0, dropping linearly to 
50% of the sky at £;sfd - 0.2. 

Two effects determine the dependence of typical color 
of background galaxies on ^sfd for a given reddening 
law. Primarily, the population of galaxies will redden (or 
'bluen') as £^sfd increases, dependent on the differences 
in extinction between the bands in question. There is a 
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Figure 3. Reddening measurement of the G ALEX- WISE galax- 
ies as a function of £^SFD- The top panel shows the median color of 
galaxies as a function of ^sfd, as well as the population reddening 
effect of Equation O The bottom panel shows a second-order poly- 
nomial fit to the reddening corrected for the population reddening 
effect (see Equation [6]), representing the expected reddening of an 
object (black) and its derivative, the value of Rfuv —Rnuv (red). 

secondary effect that stems from the extinction of galax- 
ies out of the sample population. If a population does 
not have a fixed color as a function of the magnitude on 
which it is selected, the typical color of the galaxies will 
change simply by removing the dimmest galaxies from 
the sample in regions of higher ^sfd , which we call the 
population reddening effect. This is to say 



{a-h) = {a- b)o -^[{Ra- Rb) + R. 



d{a - b) 



Es 



FD, 



(3) 

where {a — b) is the observed color of a reddened galaxy 
population in a — b color, (a — b)o is the unreddened 
color, Ra is the extinction in band a per unit £^sfd, 
and c is the extinguished band upon which the sample 
is cut at Cmax- Often, Ra — Rb is of the order of Rc^ 
and ^^^^ << 1. This is the case for the PGIO galax- 
ies, which are measured in the optical and standardized 
to have fixed color with magnitude, and thus the pop- 
ulation reddening effect is negligible. In the G ALEX- 
WISE sample, Ra - Rb = Rfuv - Rnuv ^ -0.02, 
Re = Rnuv — 8, and the population is uncorrected as a 
function of magnitude. Thus, the population reddening 

effect is non- negligible. In practice, we measure r}^~ ^ 
by selecting a large number of galaxies from our sam- 
ple with extremely low extinction, artificially extinguish 
them, re-apply the selection criterion, and measure the 
change in the median color of the population. We em- 



pl oy he re the same comparison sample (CS) as described 
ingp21( 0.02<£;sFD <0.03). 

Figure [2] shows the results of the extinction analysis. 
The £^SFD cumulative distribution function (CDF) in the 
top panel demonstrates that half of the entire sky has an 
extinction below £^sfd = 0.07, showing the importance 
of extinction parameters at these low values of £^sfd • We 
also show the CDF of the G ALEX- WISE galaxies against 
^SFD, making plain the very small fraction of our 373,303 
galaxies we observed towards £^sfd > 0.2. The second 
panel shows the number density of galaxies in the sam- 
ple, clearly showing the effects of extinction, along with 
the number density of the CS. The anomalous turn-over 
below £^SFD = 0.02 is discussed in the Appendix. The 
third panel shows Rnuv as a function of ^sfd- These 
values are measured by comparing the number density 
of galaxies in the CS to the values within a given bin. A 
second-order polynomial fit is shown. 



Rnuv = (8.36 ± 0.42) + 

(14.3 ± 6.9) ^sFD + (-82.8 ± 26) ^|fd, 



(4) 



along with the 95% formal confidence interval of the fit. 
Figure [3] shows the results of the reddening analysis. 
The first panel shows the median color of the galaxies 
in each bin, the left hand side of Equation [3l Errors 
are determined as a maximum of bootstrap and Pois- 
son errors as in ^3.11 Also shown is the population red- 
dening effect we discussed above, {FUV — NUV) cs + 
Rnuv onuv ^^sfd, as measured in the CS. We ex- 
trapolate this trend down to £^sfd = 0. The second panel 
shows a second-order polynomial fit to 

{FUV - NUV) - {FUV - NUV)o 
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{Rfuv — Rnuv) Fsfb = 

(2.11 ± 0.08) £;sFD + (-5.71 ± 0.53) Ei^j^, 



(5) 



which is equivalent to the expected color change for an in- 
dividual object due to reddening. We can determine the 
difference between the reddening coefficients, as shown 
in the bottom panel in red: 



Rfuv — Rnuv = 

(2.11 ± 0.08) + (-5.71 ± 0.53) £;sfd. 



(6) 



Both fits are shown with 95% confidence intervals. Also 
shown are the expectations from 0'D94 and F99 extinc- 
tion laws. Taking the sum of Equation S] and Equation 
[6] we find 



Rfuv = (10.47 ± 0.43) -f 

(8.59 ± 6.9) £;sFD + (-82.8 ± 26) £;|fd. 



(7) 



As a sanity check we apply the same analysis to the 
sample of 95,846 QSOs selected in ^2.31 These results are 
in agreement with the reddening analysis of the G ALEX- 
WISE galaxy sample, reproducing both the shape and 
amplitude of the curves shown in the bottom panel of 
Figure [3l albeit with larger scatter. 

Neither our observed extinction nor reddening are con- 
sistent with the literature. The standard values for ex- 
tinction are Rnuv = 8.3,8.1 and Rfuv — Rnuv = 
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Figure 4. Maps of £^sfd (top) and galaxy color (bottom). The right side columns show the PGIO galaxies, the left side columns show 
the GALEX-WISE galaxies. Maps are shown as a pair of hemispheres, Galactic North (left) and Galactic South (right). Maps are limited 
to the area in which the noise in the optical and UV color maps is less than 0.005 and 0.045 magnitudes, respectively. 



-0.026,-0.022 for the typical Ry = 3.1 in the 0'D94 
and F99 extinction laws, respectively. Our values for 
NUV extinction are inconsistent for the vast majority of 
our sample, and our FUV — NUV reddening values are 
inconsistent throughout. We note that neither of these 
values is consistent with the recalibration of SFD pro- 
posed by SFll, who predict a value of Rnuv = 6.86. 
SFll note that for £^sfd < 0.2, the SFD normalization 
is detect ably higher, closer to the original value proposed 
in SFD, but do not find any evidence for a normalization 
greater than unity that would be needed to produce the 
offset we see in Rnuv- 

4.2. Color Maps 

Figure H] shows the median color Cpuv-Nuv^ and 
£^SFD maps generated the method described in §3.11 for 
the UV galaxies. Cfuv-nuv varies dramatically over 
the map, from 0.35 to > 0.8. There does not seem to be 
any clear correlation between Cpuv-Nuv and £^sfd ex- 
cept that redder UV colors only occur in regions of high 
extinction, whereas bluer colors can be found in regions 
of both high and low extinction. 

4.3. Contrast with Optical Reddening 

Reddening laws make predictions for the extinction as 
a function of wavelength as a function of some measure 
of total dust column (here we use E'sfd), and some other 
parameter that characterizes the grain distribution, typ- 
ically Ry ' To compare to reddening laws from the litera- 
ture, we restrict our examination to the area of sky where 
we have both UV and optical extinction measures, the 
SDSS area. We show the distribution of £^sfd against 
galaxy extinction for both UV and optical in Figure [H 
In the left panel we show the expected reddening for 



these wavebands, using both 0'D94 (nearly equivalent 
to CCM in these filters) and F99 reddening laws and 
a range of Ry- It is clearly evident that the canonical 
Ry = 3.1, broadly consistent with the optical data using 
either reddening law, cannot explain the distribution of 
galaxy colors in the UV. The UV over this area is best fit 
with a F99 reddening law with Ry c^ 2.2. As a further 
check we compare to re ddening laws in the Magellanic 
Clouds as measured by (jGordon et al.|[2QQ3[ ). None of 
these reddening laws fit our data either, as shown in the 
right panels of Figure [5l Finally, we compare the high 
latitude UV extinction properties to each of the redden- 
ing la ws toward 417 stars as measured by (|Valencic et al.l 
l2004l ). Of these 417 stars, only two marginally match the 
measured Rpuv and Rnuv as measured at a typical high 
latitude value of £^sfd = 0.05 (Equation [Hand Equation 
[6]). These stars, HD 13659 and HD 204827, have mea- 
sured Ry = 2.53, 2.58, respectively, once again incon- 
sistent with well-tested value of Ry = 3.1. Thus, high 
latitude extinction is qualitatively inconsistent with typi- 
cal reddening laws in the Galaxy, the Magellanic Clouds, 
as well as sighthnes to many individual Galactic stars. 

It was shown in PGIO that the deviation in Cg-r from 
the expected value was significant in many regions of the 
SDSS footprint. To investigate the relationship between 
this variation and the equivalent variation in the UV, we 
show in Figure [6] the residual C for the GALEX-WISE 
and PGIO galaxies after subtracting the median trend 
with £^SFD shown in Figure [5l We find that the devi- 
ation in the UV color, /S.Cfuv-nuVi is also in excess 
of what would be predicted from the errors, and that 
there is no detectable correlation between /S.Cfuv-nuv 
and ACg_r. This indicates that whatever ISM process is 
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Figure 5. £^sfd versus galaxy color, C . Top left panel shows data from the G ALEX- WISE galaxy sample, bottom left shows data 
from SDSS PGIO galaxy sample, both restricted to the SDSS footprint, for comparison. The gray lines show the median galaxy color as a 
function of ^SFD- Also shown are 0'D94 (dotted) and F99 (dashed) reddening curves with Ry = 2.2, 3.1, and 5.5 (blue, green, and red, 
respectively). Neither reddening law is capable of reproducing t he color variation i n both the optical and UV. On the right the same data 
are shown with four Magellanic Cloud extinction curves from ([Gordon et al. I [20031 ). Similarly, none of these extinction curves match the 
data. 



driving variation in Rfuv — Rnuv across the sky is not 
the same as the process driving Rg — R^. Specificahy, it is 
not simply errors in ^sfd generating both effects — the 
reddening curve at high latitude is both different than 
predicted by 0'D94 and F99, and its variations also do 
not follow a simple one-parameter model. 

4.4. Extinction Modeling 

We have shown that reddening in the UV over the 
GALEX sky, Rfuv — Rnuv^ does not follow the canon- 
ical laws in 0'D94 or F99, and that variation in Rfuv — 
Rnuv on the sky is neither simply a function of error 
in E'sFD nor driven by whatever grain population vari- 
ation that generates the variation in the optical. To 
investigate the variation of UV extinction across the 
sky we model the full observed GALEX- WI SE g alaxy 
distribution using the method described in ^3.21 We 
parameterize each model by ARjjy = Rfuv — Rnuv 
and Rjjy = {Rfuv + Rnuv) /2, as these quantities 
are largely independent; the first is largely driven by 
Cfuv-nuv and the second by N^uv and Nfuv- Af- 



ter exploring the wide range of models we find we can 
fit the bulk of the data points with a grid of 15 mod- 
els, five values of AR\jy = [—1.5,0, 1.5,3,4.5] and three 
values of ^uv = [7,9, 11]. The results of the modeling 
are shown in Figure [3 We find that there are coherent 
ARjjy structures on the sky that span radians in extent. 
The northern hemisphere seems to have more regions of 
large Ai^uv, while the southern sky has more regions 
of low Ai^uv, with a number of regions showing consis- 
tent "bluening". At high latitude, where £^sfd ^ 0.04, 
we cannot discernbetween Ai^uv models. We have less 
success modeling Rjjy^ as we are very susceptible to vari- 
ation in galaxy counts due to LSS. We do find large re- 
gions of jthe sky somewhat more prone to lower R\jy and 
higher i?uv, and these regions do not seem to have a 
direct correspondence to the variations in Ai^uv- 



5. DISCUSSION 
5.1. Implications for UV measurements 
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Figure 6. The deviation of Cg_r vs. the deviation oiCpuv-NUV 
over the SDSS DR7 footprint. Deviations are measured from the 
roUing median shown in gray in each panel of Figure \E\ The red 
eUipse encloses 50% of the data points, the blue ellipse would en- 
close 50% of the data points were they all to be distributed from 
zero only according to their errors (i. e. have a reduced x^ = 1 
to a model in which the data were simply scattered about the me- 
dian according to their errors). The median error bar is shown 
in the upper left. We see that the scatter is significantly larger 
than would be expected from the error bars, and that there is no 
detectable correlation between ACg_r and ACpuy^Nuv This 
indicates that variation in reddening from the typical value in the 
UV is not primarily driven by the same ISM process that creates 
the variation in the optical reddening. 

Given that all UV extragalactic studies require a 
Galactic reddening correction, we might ask what a cor- 
rection to the standard prescription will have on these 
studies. While many analyses will be insensitive to the 
small corrections we have measured, dust attenuation 
measures using FUV — NUV colors (e.g., Calzetti et al. 
[2000), and counts analyses (e.g. number counts, corre- 
lation functions) will often be quite sensitive to Galactic 
corrections on the scale we have detected. 

Dust attenuation measuren ients using the UV SED 
slope follow the spirit of the iMeurer et al.l (|1999f ) pre- 
scription, and in general find an approximate depen- 
dence Apuv ^ 4:{FUV — NUV) with considerable 
scatter in the average slope and in the intrinsic dis- 
tributi on of the relations hip across the galaxy popu- 
lation (jSeibert et al.l l2QQ5f ) — despite the many caveats, 
this measure has been applied to UV luminous galax- 
ies across a wide range of redshifts. Our results imply 
a correction to -A(FUV-NUV) -0.05-0.1 (for a typi- 
cal extragalactic line of sight with E(B-V) - 0.03-0.07) 
implying a systematic dust attenuation correction of 
—AAFUY^cor ^ 0.2 — 0.4 (with galaxies intrinsically bluer 
than previously determined). The effect on counts and 
correlation function analyses will depend on the slope of 
the extragalactic counts and fiuctuations in the Galactic 
dust column/extinction correction. The impact of Galac 



tic reddening has been discu ssed pre viousl y (IXu et al. 
|2005; Milhard et al. 2007; Hei nis et al.ll2009l : lBlake et al. 
[2OIO1 : iHammer et al.l i2010l) . Here we merely emphasize 
that the implied systematic correction will provide an ad- 



ditional source of hitherto uncorrected variance for the 
counts analyses and component of the measured power 
spectrum that remains highly correlated with Galactic 
structure. 

Galactic reddening corrections may also impact QSO 
colors, which may then be interpreted as a variation in 
intrinsic QSO reddening, or absorption of UV photons 
by intergalactic HI or dust. As discussed above, we 
have tested our methodology using a smaller QSO sam- 
ple and find similar results as with our main galaxy sam- 
ple. Observed- frame UV measurements have been used 
to determine the rest-frame EUV/UV spectrum of QSO s 
(jZheng et al.l [19971 : Piblfer et al.l[200l iShuh et al.l 120121 ). 
requiring accurate Galactic reddening corrections. Pre- 
cise measurements may be able to detect the "Lyman 
valley", Lyman continuum absorption resulting from 
the integra ted effect of many low density absorbers. 
iMoller fc Jakobsen (1990i) predict Lyman valley absorp- 
tion at the rsj 5% level (at 700A) out to 2: ~ 1, an ampli- 
tude comparable to the Galactic reddening corrections 
proposed here. In the most com prehensive study o f the 
EUV spectrum of QSOs to date, iShull et al.l (j2012r ) are 
unable to detect this valley, though just at the limits of 
their statistics and using standard Galactic corrections. 
Improved reddening corrections will better constrain this 
valley and the slopes of the ultraviolet composite spec- 
trum in the EUV. 

5.2. Implications for Galactic Dust 

Reddening in the UV is thought to be governed by two 
groups of particles; molecule sized carbonaceous grains 
including PAHs and v e ry sma ll silicate grain s (VSGs; 
iWeingartner fc Draind l2001bl . iDesert et all I1990L 
Zubko et al. 2004). While the details of the models vary, 
PAHs and associated grains typically have a stronger ex- 
tinction in the NUV (near the 2175 A bump) and VSGs 
have more extinction toward the FUV. Some models re- 
fiect an additional large grain component responsible for 
reddening in the optical and relatively gray extinction 
in the UV, while others lump all silicate grains together 
with a grain size distribution. As standard extinction 
laws of the Milky Way parameterized only by Ry fail 
in the UV for the low extinction areas surveyed in this 
work (Figure [5]), we are interested in which components 
of the physical dust models could be altered to match our 
observations. Specifically, we wish to match the results 
that 1) at high latitude extinction in the NUV is ^ 10% 
larger than expected, extinction in the FUV is ~ 30% 
larger than expected, and reddening in the optical is un- 
changed. 

Modifying the strength of any large grain component 
will do relatively little to match our observations. While 
we do see a 10% increase in Rnuv at £^sfd < 0.15, we 
see no change in optical reddening. Since larger grains 
strongly affect reddening in the optical, variation in the 
larger grains is unlikely to be the dominant effect. The 
carbonaceous-PAH component typically provides more 
extinction in the NUV than the FUV GALEX bands, 
thus variation in this component provides correlation be- 
tween increased overall extinction and decreased UV red- 
dening. As this is the opposite of the effect we see, varia- 
tion in the carbonaceous-PAH is unlikely to be the major 
culprit. Varying the VSG component may be the best 
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Figure 7. Mapsof the goodness-of-fit of each of our 15 Monte Carlo models. Top panels represent ARjjy = Rpuv — Rnuv ^ bottom 
panels represent Rjjy = (Rfuv + Rnuv) /2; Left panels show the northern galactic hemisphere, right panels show the southern Galactic 

hemisphere. The color of each 16 square degree region is colored according to the goodness-of-fit statistic, yXred^ ^^ ^^^ ^^^ each model. In 
the top panels, red, green, and blue correspond to Ai^uv = [0? 1-5, 3], with the brightness in each color representing the 
fitting model across the range of ^uv- In the bottom panels, red, green, and blue correspond to ^uv = [7, 9, 11], with the brightness in each 
color representing the yXred ^^ ^^^ ^^^^ fitting model across the range of ARuy. Additionally, in the top panels, if the ARjjy = [0, 1.5, 3] 

models fail, each having a best-fit yX^ed ^ ^' ^^^ ^^^^ ^^ among ARjjy = [—1.5,4.5] is shown as a circle in blue or red respectively, with 
intensity of color once again representing goodness of fit. 

The maps show that there are coherent radian-scale features in ARjjy over largeswaths of the sky, with a bias toward higher ARjjy in 
the north and lower ARjjy in the south. Models are less sensitive to variation in Rjjy, but we do detect some large structures. Structures 
in the two metrics do not seem to be correlated. 

option to fit our observations, providing a positive corre- 
lation between reddening and extinction in the UV. We 
note that since the extinction dramaticahy rises from the 
NUV to the FUV (10% to 35%) as compared to stan- 
dard IVIilky Way reddening laws, and the optical redden- 
ing is unaffected, there must be a relative variation in 
the smallest silicate grains compared to larger silicate 
grains, rather than an overall renormalization of the sili- 
cates. We conclude that there does not exist a universal 
grain size distribution for the IMilky Way, even allowing 
for parameterization with Ry. 
It is perhaps not surprising that the lower density en- 



vironments probed by these sightlines would harbor a 
different population of grains. Dust grains are known to 
grow in the ISM, and this growth is strong function of gas 
density (|Draineiri99Q[ ). Furthermore, grain are shattered 
in the ISIVI in shocks, including grain-grain and grain- 
ion collisions (Tielens et al. 1994). This is to say that if 
grains have been outside of the molecular environment 
long enough to be processed by shocks, we would indeed 
expect there to be an excess of small grains as observed. 
This is also in quah tative agreement with the work of 
(Par adiset aIll2Qllf ) in the LMC. We note that grain- 
size variation in the IVIW beyond the Ry prescription 
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has been noted before, in precision analyses of individ- 
ual sightlines toward denser regions (jLarson et alJ I2QQQL 
(|Clavtonet alJl2QQ3[ )). 

We also would expect there to be significant variation 
in the the fraction of PAHs as we look towa rd low density 
envir onments. Obse rvations of the LMC (Paradis et all 



I2QQ9D and the SMC (jSandstrom et alJl2QlQ[ ) both show a 
correlation between the fraction of dust in PAHs and the 
density of the environment bei ng probed. In our own 
Galaxy (jGillmon fc Shulll l2QQ6f ) showed that there is a 
significant decrease in molecular fraction at column den- 
sities represented by £^sfd < 0.1. While VSG variation 
seems to dominate PAH variation in aggregate, we do 
see evidence for a variation in Figure [7] that cannot be 
explained by only varying the VSG fraction. Were all 
variation to be limited to the VSG fraction, there would 
be a strong positive correlation between UV color (top 
of Figure [7j) and UV extinction (bottom of Figure [7j) . As 
this is not the case, there must be some added variation 
in dust grain properties at intermediate latitudes with an 
anti-correlation between reddening and extinction. This 
effect is consistent with variations of PAH fraction across 
the high latitude sky. 

By connecting the angular structure in the ARjjy 
maps (Figure [7j) to the angular structure of other tracers 
of ISM structure, we could get a better physical under- 
standing of the drivers of grain shattering and PAH frac- 
tion variation. A brief visual inspec tion of the location s 
of the four high-latitude radio loops (Hasla m et al.l ll97?). 
the H-alpha sky (F inkbeiner 2003) . the soft X-ray back- 
ground skv (iFrevb erg 1998 1), and the HI sky as a function 
of velocity (jKalberla et al.ll2005f ) shows no convincing as- 
sociation with the ARjjy maps. It is not entirely clear 
what to make of this lack of association, other than it 
may help to rule out some very simple methods of grain 
destruction that might predict a correlation. A more de- 
tailed comparison would be valuable, but is beyond the 
scope of this work. 

6. CONCLUSIONS 

Simply put, ultraviolet extinction at high latitude is 
neither quantitatively nor qualitatively consistent with 
the standard Galactic reddening laws. More specifically: 

1. The overall ultraviolet extinction parameter 
Rnuv — 9 for £^SFD — 0.1, trending toward 
Rnuv — 8 near £^sfd — 0.2 (see Equation H]). The 
UV reddening can be parameterized as in Equation 
[HI starting at Rpuv — Rnuv ^ 2 and dropping to- 
ward as £^SFD increases. At low extinction this 
represents Rpuv and Rnuv at 1.3 and 1.1 times 
their fiducial values and 1.6 and 1.3 times the val- 
ues suggested by the rescaling of SFD in SFll. 

2. The premise that dust extinction across wave- 
lengths can be parameterized along the line sight 
by a single parameter, i?y, is not consistent with 
observations at high latitude. We find that over 
the SDSS spectroscopic area optical data are con- 
sistent with Ry = 3.1. In the same region, the 
UV reddening is inconsistent with CCM or 0'D94 
reddening, and only consisting with F99 redden- 
ing for Ry c^ 2.2. Furthermore, both optical and 
UV reddening depart from these Ry values signifi- 



cantly over the observed area, but these variations 
are not correlated with one another. 

3. The UV high latitude sky is neither consistent with 
a single reddening parameter, nor a range of red- 
dening parameters characterized by £^sfd- Instead, 
AR\jY = Rpuv — Rnuv ranges from -1.5 to 4.5, 
in radian-scale coherent structures across the sky. 
Variation in R\jy is less clear, but there seem to 
be large, coherent structures that do not follow the 
structures in ARjjy. 

4. Extragalactic studies that use FUV — NUV colors 
as dust attenuation measures, counts and correla- 
tion analyses, QSO EUV spectral energy distribu- 
tion and low redshift {z < 1) IGM absorption mea- 
surements may be sufficiently sensitive to require 
consideration of the variation in Galactic extinc- 
tion and reddening. 

5. These results are consistent with low density, high 
latitude ISM having a typically higher fraction 
of very small silicate grains, which generate both 
higher extinction and stronger reddening in the 
GALEX UV channels. 
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Appendix 

Our analysis presupposes that the positions of galaxies 
on the sky are not correlated with the foreground extinc- 
tion in the Galaxy. While this is likely true in reality, 
throughout this work we use the SFD extinction map 
to determ ine foreground extinction, which may itself be 
incorrect. lYahata et al.l ()2007l :Y07) purport to find just 
such a correlation which we must treat with great care. 
They showed that the number counts of galaxies found 
photometrically in SDSS DR4 were significantly lower 
than expected at very low £^sfd- They develop a conjec- 
ture that the contribution of far-IR fiux from the galax- 
ies themselves create this effect by artificially increasing 
£^SFD along lines of sight to galaxies, thus decreasing the 
number of galaxies towards low extinction. The fact that 
galaxies contribute significantly to the SFD map is con- 
firmed in .Kashiwagi et al.. (^20 12.) . 
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To test this effect in our data, we compare the £^sfd to- 
ward our galaxies to the Esfd for positions offset by ±A/. 
If the galaxies do indeed contribute significantly to £^sfd , 
there should be a detectable increase in £^sfd towards the 
galaxies as compared to all the positions slightly offset 
in Galactic longitude. In Figure [8] we see just such an 
effect in both the PGIO and GALEX-WISE populations. 
The difference rises steeply within the ~ 6' ^sfd pixel, 
and flattens by about A/ = 25^ The effect is weaker for 
the PGIO galaxies than the GALEX-WISE galaxies, and 
increases in strength when only a subset of the brighter 
GALEX-WISE galaxies is selected. All of these assess- 
ments are qualitatively in agreement with the conjecture 
in YQ7 and quantitatively in agreement with the results 
of iKashiwagi et ahl ()2Q12f ). o To determine whether this 
effect is strong enough to produce the distribution of 
galaxies we observe, we build a simple model. We take 
the CS galaxies and artificially extinguish them, remove 
galaxies with NUV > 20, and bin them by extinction. 
As shown in Figure [H this model roughly reproduces the 
exponential slope in galaxy number density vs. ^sfd for 
Rnuv = 8.5. The model does not reproduce the anoma- 
lous turn-over at £^sfd = 0.02 in our data, which is anal- 
ogous to the effect observed first in Y07. We also show 
the same model, but now include the detected correla- 
tion between GALEX-WISE galaxies and excess £^sfd as 
a function of galaxy NUV magnitude. While there is a 
marginal effect on the model at extremely low extinction, 
it is far too weak to reproduce the turnover. We conclude 
that the decrement of galaxies at low extinction, while 
clearly visible in our data, is not due to galaxy-induced 
variation in £^sfd- We do not model any overall offset in 
the £^SFD zero-point that might stem from the aggregate 
of all galaxies. While such an effect could reduce the sig- 
nificance of the SFD signal toward low £^sfd, and thus 
flatten the slope of galaxy number density histogram, we 
do not find any way that this effect can create a turn- 
over. It remains unclear what causes the decrement in 
galaxies where £^sfd < 0.02. To avoid contaminating our 
analysis with this clear error, we ignore the few areas of 
sky where £^sfd < 0.02 when galaxy number density is 
under consideration. 

Y07 is concerned with systematically incorrect fore- 
ground extinctions toward galaxies, and the impact of 
that error on the number density of background galax- 
ies. Far more nefarious for our data is a known intrin- 
sic dependence of ga laxy UV color on far IR luminosity 
(|Calzetti et al.l 120001 ). Galaxies that are redder in the 
UV tend to have stronger far IR flux, which could bias 
our results if that same far IR light systematically con- 
taminates £^SFD- To address this effect we investigate 
the variation of median UV galaxy color as a function 
of the this erroneous IR component associated with the 
galaxies. We examine the distribution of A£^sfd [25^], 
the typical increase in £^sfd at a galaxy as compared 
to positions offset by ±25' in Galactic longitude. The 
distribution is shown in Figure [3l We then examine 
median galaxy color, C, as a function of A£^sfd [25^]. 
We have found that galaxies typically become somewhat 
redder in the UV as £^sfd increases (see Figure [3]), and 
since variability and amplitude in £^sfd are correlated, 
there is a spurious correlation between galaxy color and 
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Figure 8. The change in median value of ^SFD as we dis- 
place Galactic latitude away from our galaxies, ^sfd [I] — 
(^SFD [I + ^l] + ^SFD [I — ^l\) /2. We do detect contamination 
at the sub-millimagnitudesub-millimag level in both PGIO and 
GALEX-WISE galaxies. This effect is consistent with far-IR as- 
sociated with the galaxies themselves, though we find it to be too 
small to accou nt for the observed d ecrement of low -Esfd galaxies 
as proposed in lYahata et al.l (|2007l '). 

|A£^SFD [25^] I . To nullify this effect we examine the dif- 
ference in median color of galaxies between excess and 
decrement IR, 



^Cfuv-nuv (A£^sfd [25']) = 
Cfuv-nuv (A£^sfd [25 ]) — 
Cfuv-nuv (— A£^sfd [25 ]) . 



(8) 



In the bottom panel of figure [10] we show that there is a 
trend in galaxy UV color with A£^sfd [25^] of ~ 2 magni- 
tudes per magnitude variation in A£^sfd- We find that 
a distribution of galaxies with this observed A£^sfd and 
UV color trend only generates a typical shift in UV color 
of about 0.6 millimagnitudes, far below our sensitivity. 
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